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Abstract. Thermophoresis is a tool often applied in complex plasma 
experiments. One of the usual stated benefits over other experimental tools is 
that changes induced by thermophoresis neither directly depend on, nor directly 
influence, the plasma parameters. From electronic data, plasma emission profiles 
in the sheath, and Langmuir probe data in the plasma bulk, we conclude that this 
assumption does not hold. An important effect on the levitation of dust particles 
in argon plasma is observed as well. The reason behind the changes in plasma 
parameters seems to be the change in neutral atom density accompanying the 
increased gas temperature while running at constant pressure. 
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1. Thermophoresis in complex plasma 

Partially ionized gases containing small solid particles, called complex plasmas, are a 
widespread phenomenon in space, found for instance in molecular clouds, planetary 
atmospheres, planetary rings, and comets [2 [2]- In the laboratory, complex plasma 
has become a promising tool to study the physics of solid state systems, fluid dynamics 
and turbulence, even biophysics, on a level accessible with ordinary optical techniques 
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Unless experiments are performed during parabolic flights, or under micro- 
gravity conditions on the International Space Station [6], gravity dominates the force 
equilibrium in complex plasma, resulting in the formation of two-dimensional crystal 
structures in places where the electrostatic force acting on the particles balances 
gravity. An alternative is to extend the confinement of the particles vertically 
by adding a dielectric (glass) box in the discharge [7]. In these experiments, 
thermophoresis played an important role in reducing the effect of gravity, and 3- 
dimensional crystals were produced. 

The thermophorctic force depends on the establishment of a temperature gradient 
in the background gas, resulting in a net imbalance of the momentum imparted by 
neutral atoms colliding with a particle in the gas. Therefore, it has to date been 
regarded as a tool independent of the plasma parameters, which in complex plasma 
experiments are hard to measure, since they can become strongly coupled to the 
presence of particles in the plasma [8]. As a consequence, changes in discharge 
characteristics (like the DC bias, or absorbed power) accompanying temperature 
changes are usually considered insignificant and are not specified in complex plasma 
experiments involving thermophoresis, see for instance [9]. 

This paper reports on complex plasma experiments involving thermophoresis 
performed in a modified Gaseous Electronic Conference (GEC) reference cell [TU] at 
the Center for Astrophyiscs, Space Phyiscs, and Engineering Research (CASPER). 
Two direct observations connected with thermophoresis are a change in the DC bias 
on the powered electrode with increasing electrode temperature, which was briefly 
mentioned before [11], accompanied by a change in levitation height of the particles. 
From electronic data, optical emission data, and Langmuir probe data, we show that 
the plasma properties indeed change. Using a self-consistent fluid model [12) . we 
show that the volume averaged gas-temperature increases with electrode temperature, 
resulting in a decrease of the neutral atom density. 

2. The complex plasma experiment at CASPER 

The experiments discussed here were performed using argon gas at low pressures (a few 
tens of Pascals). The upper electrode in the cell, as well as the walls surrounding the 
chamber, are grounded, while a 13.56 MHz radio-frequency (RF) potential is applied 
to the lower electrode with an amplitude of a few tens of Volts. 

Once plasma is formed, spherical micron-sized melamine-formaldehyde (MF) 
particles are introduced from shakers situated at the top of the chamber. The particles 
fall through the grounded electrode, which consists of a hollow cylinder, rather than 
a solid electrode. At equilibrium, the particles are vertically suspended in the sheath 
above the bottom electrode, while horizontally confined by a cutout plate containing 
a cylindrical depression creating a radial electric field placed on top of the bottom 
electrode. 
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Two diode lasers are used to illuminate the particles. The lasers are equipped 
with cylindrical lenses, so that thin (±100 micron), but wide (fan angle > 5 °) laser 
sheets are created. A vertical laser sheet in combination with a side-view camera 
allows side-view pictures of the suspended particles, whereas a horizontal laser sheet 
together with a top- view camera focused through the hollow upper cylinder, allows for 
top- view images. A schematic drawing of the experimental setup is shown in figure [T] 



The lower electrode temperature is regulated by coolant liquid running through 
tubes which are placed inside of bores in the electrode. The temperature of the liquid is 
maintained by a heater/chiller with a thermostat, and is adjusted from 0° C to +70° 
in steps of 10° C. The lower electrode temperature is measured with a thermistor 
attached to the bottom of the electrode and logged with a digital datalogger. 

Electronic data, including the driving potential, the electrode current, and the 
signals from two derivative probes, measuring the derivatives of the potential and 
electrode current, are collected on two oscilloscopes. The latter signals are used to 
determine the phase angle between the potential and the current. From this angle 
and the first two signals, the power input into the plasma is calculated. From the 
electronic data, the floating DC bias on the lower electrode is obtained too, as the 
shift in the mean value of the driving potential data. 

3. Data obtained during thermophoresis complex plasma experiments 

This section presents the data obtained during experiments at a pressure of 200 mTorr 
(26 Pa). In these experiments 8.89 micron diameter MF particles were used. Initially, 
the DC bias was allowed to float. Upon observing the change in DC bias with electrode 
temperature, the experiments were repeated using an external DC power supply to 
fix the DC bias. We first present the measured DC bias and the levitation height of 
the particles above the lower electrode obtained with the side-view camera, for the 
two sets of experiments. We then present the additional data obtained from the data- 
loggers, oscilloscopes, and Langmuir probe, as well as optical emission profiles in the 
sheath. 





CCD Camera 
Vertical Sma^e 



Figure 1. A sketch of the experimental setup used. Shown are the laser/camera 
pairs, the lower electrode with the cutout milled in the cover plate and a cloud of 
dust particles suspended above the cutout. 
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3.1. DC bias and dust levitation 

The measured DC bias on the lower electrode is shown for both sets of data in figure 
[21 When the DC bias is allowed to float, it becomes less negative with increasing 
lower electrode temperature. Because our discharge chamber is fairly symmetric, the 
induced DC bias is small; even so, the change in DC bias with the change in lower 
electrode temperature is significant. Even when the power supply was used, the DC 
bias had to be manually adjusted for each electrode temperature setting, showing that 
it still changed despite the connection to the power supply. Overall, this resulted in a 
slight variation of the DC bias, although the variation was much smaller than in the 
experiments with floating DC bias. 
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Figure 2. The mean floating DC bias measured for thermophoresis experiments 
with different thermostat settings. The bars indicate the range of bias potential 
measured during the 150 second measurement period, for each lower electrode 
temperature setting. X refers to the floating DC bias experiment, while the 
refer to the experiment with the DC bias held constant by the external DC power 
supply. 



For each temperature setting in each experiment, the DC bias was logged for 150 
seconds. The error bars in figure [2] indicate the range over which the DC bias varied 
during these 150 seconds, whereas the symbol shows the mean. As can be seen, the 
DC bias varies over a rather large range. This is better illustrated in figure [3j where 
fifty consecutive data points are shown. The lower panel shows the registered lower 
electrode temperature with the thermostat set at 70° C, while the upper panel shows 
the DC bias on the lower electrode. In this case, the DC bias was floating with values 
fluctuating around -5.45 V. Clearly, the DC bias directly follows the temperature on 
the powered electrode. The variations in temperature shown are due to hysteresis 
in the heater/thermostat system. As can be seen, there is an offset between the 
temperature setting on the thermostat and the mean temperature measured with 
the thermistor (about 4° C at 70° C), as well as the tendency to overshoot the set 
temperature (plus offset). 

The levitation height of the particles obtained in both sets of experiments is shown 
in figure H]B, while in figure A the overlayed side-view images for both experiments 
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Figure 3. The electrode temperature and DC bias measured over approximately 
150 seconds for one of the thermophoresis experiments described in the text. 
The DC bias clearly changes with electrode temperature over shorter time-scales, 
illustrating that the change in DC bias is not a long-term charging effect. 



at an electrode temperature of 40° is shown. The white dots represent the particles 
in the experiment with fixed DC bias, the black dots the particles in the experiment 
with floating DC bias. There is a clear difference in levitation height of about 20 
pixels. The sideview camera resolution is roughly 25 (im per pixel, so this difference 
in levitation height corresponds to 500 fim. 
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Figure 4. A. This image shows the result when a side-view image of dust particles 
levitated with floating DC bias was subtracted from an image of dust particles 
levitated with fixed DC bias, with the electrode held at 40° C. The black dots 
correspond to the dust particles with floating DC bias, the bright dots with fixed 
DC bias. Clearly, a difference in levitation height is observed. B. The mean 
levitation height for experiments with fixed DC bias (X) and the experiments 
with floating DC bias (o). The size of the symbols is an indication of the spread 
in the mean levitation height obtained from tens of images taken during each 
experiment. 



At the lowest electrode temperature setting (0° C) the dust levitation height 
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is the same for both sets of experiments. For increasing temperature the difference 
increases, coinciding with the increasing difference in the DC bias between the two 
sets of experiments. At the highest electrode temperature setting (70 °), the difference 
in levitation height decreases again. 

3.2. Electronic data from oscilloscopes 

To investigate whether the changes shown are related to changes in the plasma 
parameters, we logged the root- mean-square driving potential (Vrms) and current 
(Irms), as well as signals from two derivative probes. The raw signals are shown 
in figure [5] Each RF cycle is sampled by 37 points, giving a resolution of roughly 2 
nanoseconds. The top panel shows Vrms and 300 x Irms (which is typically on the 
order of 10 mA to fOO mA). The bottom panel shows the normalized derivative probe 
signals, which are used to determine the phase angle $ between Vrms and Irms- 




Data point (x 2 microseconds) 
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Figure 5. Electronic data obtained from the oscilloscopes. The top panel shows 
the RMS potential on the lower electrode, as well as the net current to the lower 
electrode (x 300). The bottom panel shows the normalized derivative probe 
signals, which are used to obtain the phase angle $ between the current and the 
potential. 

Using this data, we can calculate the total power absorbed in the plasma, as 
Prms = VrmsIrms cos($). The amplitude of the RMS voltage and current, the 
phase angle as well as the RMS plasma power derived for both sets of experiments 
are shown in figure [6l From the results, we see that the amplitude of the RMS 
potential decreases, the amplitude of the RMS current increases and $ decreases, 
as the electrode temperature is increased. The final result is an increase in the RMS 
plasma power, with a linear fit to both sets yielding an increase in power of 2 mW/°C, 
so that over the temperature range considered, the power increases by 10%. 

3.3. Bulk Langmuir probe data 

Using a SmartProbe Langmuir probe from Scientific Systems |13j the plasma 
parameters in the bulk plasma were measured while varying the electrode temperature. 
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Figure 6. The amplitude of the RMS voltage, current, the phase angle between 
them and the RMS aborbed power in the plasma, versus electrode temperature 
for the experiments with floating DC bias (□), and with fixed DC bias (•). 



Due to technical constraints with the probe, the electrode temperature setting was 
limited between 10° C and 60° C. Figure [7| shows the local floating potential (defined 
as the floating potential with respect to ground measured by the probe minus the local 
plasma potential with respect to ground as measured by the probe), the local electron 
temperature in electron volts, and the local electron density. 
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Figure 7. The Langmuir probe data obtained in the bulk plasma. From left 
to right, the local floating potential, the electron temperature and the electron 
density are shown. 



A large scatter in the data can be seen, which is often the case with Langmuir 
probe measurements. No clear trend with the electrode temperature is apparent, 
although one could argue that the electron density appears to decrease with increasing 
electrode temperature. The electron temperature does seem to be high, since in 
typical plasma discharges in the GEC cell at similar discharge settings, the electron 
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temperature is normally on the order of 1 eV [14] - 
3.4- Sheath emission profiles 

Optical emission profiles in the sheath were obtained with the side-view camera. To 
ensure that particles did not influence the profiles, they were taken in the absence of 
dust. Figure |8] shows the side-view obtained in an experiment with fixed DC bias at 
0° C substracted from the side-view obtained at 70° C. The images show the plasma 
glow fairly close to the lower electrode surface. The contrast has been enhanced in 
the image to increase visibility of the intensity differences. 

The dark region indicates that higher up in the sheath, the intensity is less at 
higher electrode temperature, while the thin bright structures show that at higher 
electrode temperature the local emission is increased in front of the electrode. These 
thin emission structures have recently been associated with the presence of fast neutral 
particles [TS] or energetic ions [16| . At the time these pictures were obtained, there 
was a slight mis-alignment between the camera and the lower electrode, which is the 
reason why the band structures are not precisely horizontal. 



Figure 8. The side-view obtained at °C substracted from the side-view at 
70° C, showing that at higher electrode temperature the intensity decreases in a 
broad region higher up in the sheath, while it increases in narrow regions right in 
front of the electrode surface. 

Figure [9] shows vertical emission profiles collected at 10, 40, and 60° C divided by 
the one obtained at 70° C. The total recorded height is 480 pixels, and the images are 
taken slightly higher up in the sheath than in figure |U The recorded hight corresponds 
to 12 mm, roughly the width of the sheath in front of the lower electrode. The profiles 
in both sets of experiments look similar, having a peak towards the bottom of the 
picture, corresponding to the highest horizontal band of figure [8] For all temperatures 
the intensity becomes the same at the top of the image, which corresponds with 
the start of the plasma bulk, indicating that the bulk plasma emission is the same 
for all temperatures. For increasing electrode temperature, the observed intensity 
throughout a large part of the sheath decreases, however, indicating a change in the 
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Figure 9. Vertical emission profiles in the sheath normalized with respect to the 
emission profile at 70 °C, for different electrode temperatures in the experiments 
with fixed DC bias (top) and with floating DC bias (bottom). 



4. Discussion of results 

4-1- Emission profiles 

Most of the visible line radiation in argon plasma at the powers and pressures 
discussed, is due to the de-excitation of neutral atoms in excited states rather than due 
to excited ionic states |17j . since the energy required for the direct excitation of ions 
is high (35 eV) and we do not expect a significant electron population with sufficient 
electron energy. The most likely excitation process is due to direct electron-impact 
excitation by energetic electrons (with E > 11.7 eV), whereas the most likely de- 
excitation process is spontaneous de-excitation with the emission of a photon. There 
are three possible explanations for the observed reduction in emission: a reduction in 
the creation of excited atoms, a different de-excitation process (which does not involve 
the emission of photons in the visible spectrum) becoming dominant, or an increase 
in transport of the excited species away from their origin. 

Due to the long life-time of meta-stable states, transport of such excited atoms 
away from their origin is possible. Since the GEC cell is constructed to minimize gas- 
flow, and since the temperature gradients induced by thermophoresis are too small to 
cause thermal convection [18], the excited neutral species will move through random 
Brownian motion. Thus, an increase in transport is possible with an increase in the 
neutral gas temperature. Since most energetic electrons are found in the (pre-)sheath, 
the main source of excited atoms is the (pre-)sheath region. Increased transport of 
these excited species might diffuse the emission coming from the (pre-)sheath region, 
explaining in part the reduced observed emission. 

Different de-population mechanisms from excited states could include the 
excitation from such a state to a higher-energy state, which in turn spontaneously 
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de-excites emitting a photon with a wavelength outside of the visual spectrum, or 
increased electron-impact ionization from the excited state. For the de-population to 
become important, a significant increase in electron density in the sheath would be 
required, since these processes require a second electron to collide with the excited 
atom within its life-time. Langmuir probe data in the sheath is not available, while 
details about the population of energy levels would require techniques such as line-ratio 
spectroscopy. We can therefore at this moment not exclude this possibility. 

Finally, the easiest explanation would sipmly be a decrease in electron density 
and/or neutral gas density, since electron- impact excitation rates are proportional 
to the product of both densities. This possibility seems more likely, as is discussed 
below. A significant cooling of the electron population as a whole would also explain 
a reduction in the excitation processes, since an electron energy of at least 11.7 eV is 
required for the most dominant excitation transition. Again, we have no direct proof 
of this. 

4-2. Langmuir probe data 

The Langmuir probe data presented above show no clear trend with temperature. 
According to numerical simulations, the local floating potential (which basically shows 
how negatively charged an object immersed in the local plasma would become) should 
be directly proportional to the electron temperature; (Vn — V p ) = —KT e where K 
is a constant that depends on the gas. At room-temperature in argon it was shown 
that K — 2.989 [19]. A plot of the ratio obtained from our data is shown in figure [TOl 
showing results consistent with theory. It also demonstrates a lack of a clear trend 
with electrode temperature. We therefore conclude that the bulk plasma parameters 
do not show a clear change with lower electrode temperature, which is consistent with 
the emission profiles. 
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Figure 10. — K determined from Langmuir probe data of the floating potential, 
plasma potential and electron temperature; —K = (Vfi — V v )/T e . The measured 
values shown are very close to the value of -2.989 as determined in numerical 
simulations 1191 . 



4-3. Electronic data 

In a low pressure RF discharge the sheaths act mostly capacitively, whereas the main 
power dissipation is due to the resistive bulk. In this case, the impedance of the 
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plasma can be represented by two sheath reactances, X s , in series with the bulk 
resistance R b [20 . Since the dissipative character of the sheaths is much smaller 
than the bulk, they are considered to be completely capacitive, and hence their 
impedance is purely imaginary. The discharge then has a complex impedance Z given 
by Z = Rb + i(X sg + X sp ) where it is assumed that the grounded sheath and the 
powered sheath are not identical. 

The absolute value of the impedance is given by \Z\ — V/I, the ratio of the 
discharge voltage and current. The phase angle between the voltage and current, $, 
is related to the resistance and reactance through tan($) = (X sg + X sp )/Rb, so that 
R b = \Z\/y/l + tan 2 ($). From our electronic data, we can thus calculate the RMS 
impedance, resistance and reactance. The results are shown in figure QT] 
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Figure 11. The magnitude of the RMS impedance Z, the resistance R, and the 
reactance X of the plasma discharge against lower electrode temperature. The 
□ corresponds to experiments with floating DC bias, the • to experiments with 
fixed DC bias. The dashed lines correspond to linear fits as noted in the graphs. 

In a capacitively coupled RF discharge in argon at the pressures and powers 
considered, the bulk resistance can be approximated by Rb = i/ c /uipeQ, with v c = 
n ga s J &v e f e dv e the electron-neutral collision frequency, a the cross section for elastic 
electron-neutral collisions [50], which is approximately 3 x 10~ 19 m 2 for the electron 
energies in the discharge [21) . and w p = \J n e e 2 / m e eo the electron plasma oscillation 
frequency. Hence, R b « 3 x I0~ 19 m e n gas y/ 8k B T e /Trm e /n e e 2 w 1CT 7 ' {n gas /n e )y/% 
Ohm, where we have assumed a Maxwellian distribution for the electrons. For a 
pressure of 200 mTorr, an electron density of about 5 x 10 14 m~ 3 , and an electron 
temperature of 7 eV, we obtain Rb ~ 370 Ohm, which is very close to the determined 
values. 

The sheath reactance can be approximated by X s = d s /uj]ipeoA sa 1321 d s /A 
Ohm, where d s is the sheath width and A is the electrode surface area [5U] • The sheath 
width can be approximated by the Child-Langmuir law, as d s « 10- 4 (AV0 3 / 4 J+ 1/2 , 
with Ay the voltage drop accros the sheath, and J+ the ion current density across 
the sheath |22j . Assuming that the ion current density can be approximated by the 
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Bohm flux (J + = en +s ^JkBT e /m + ) : and using the surface area of the lower electrode 
(A = tt(0.05) 2 ), the sheath reactance becomes X s » 10 10 (A^) 3/4 n+ 1/2 T e _1/4 Ohm. 
For a potential jump across the sheath of 30 Volts, an ion density of 5 x 10 14 m~ 3 
(assuming quasi- neutrality at the sheath edge), and an electron temperature of 7 eV, 
we find X s w 340 Ohm, which is in good agreement with the determined values. 

From the above, we see that a slight increase in resistance can result from 
an increase in gas density, a decrease in electron density, an increase in electron 
temperature, or a combination of these. On the other hand, the much stronger decrease 
in sheath reactance must result from a shrinking of the sheaths, which can result from 
an increase in plasma density, or electron temperature, or a smaller potential jump 
across the sheath. The observed change in plasma emission in the sheath might be an 
indication of these changes in the sheath width. Since the same decrease in reactance 
is also seen in the experiments with fixed DC bias, it is more likely that a change in 
the plasma parameters is responsible. 

4-4 ■ Neutral atom density 

Most complex plasma experiments, like the ones reported here, are performed at 
constant pressure. Heating or cooling the lower electrode causes heating or cooling 
of the gas inside the chamber, however. Even though one might suspect this to be a 
small effect, the change in temperature can be significant. Self-consistent fluid model 
calculations were done for our experimental setup for different electrode temperature 
settings. The volume averaged gas temperature computed with this model, which is 
used for self-consistent fluid simulations of complex plasmas |12) . is shown in figure 
[T2l Assuming the ideal gas law, P gas — n gas kBT gas , we see that the density for 200 
mTorr can vary from 6.6 x 10 21 m~ 3 at 0° C to 6.0 x 10 21 m~ 3 at 70° C, a difference 
of 10 %. Such a change in neutral atom density will have a definite effect on both 
the global particle and power balance in the discharge, in turn affecting the electron 
temperature and the electron density. 

4-5. Additional considerations 

Since the coverplate on top of the lower stainless steel electrode is made of aluminum, 
it is possible that a contact potential might be responsible for the change in DC bias. 
We therefore repeated the experiment without the coverplate and found the same 
result. This also eliminated secondary electron emission as a cause, since aluminum 
has (2-4) times the secondary electron emission coefficient of steel [23] , which means 
that any significant change in the DC bias behaviour would have been observable. 
The experiments were also repeated without the presence of particles, again showing 
the same results, indicating that the particles were not creating the observed change 
in DC bias. 

One possibility we can not exclude is a change in surface-recombination chemistry, 
due to the increased electrode temperature [24 . It is not unlikely that the surface 
accomodation might change with increased temperature, and this change might in 
part be responsible both for a change in the sheath behavior, as well as a change in 
the plasma parameters. At this time, no observational techniques are available in the 
CASPER laboratory to study this further. However, even if this were the responsible 
mechanism, it would still show that using thermophoresis in complex plasma changes 
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Figure 12. The volume averaged gas temperature in the discharge for various 
electrode temperatures, calculated using a self-consistent complex plasma fluid 
model. All other surfaces were kept at room-temperature (293 K, 20° C). The 
temperature response can be well represented by a linear function, indicated by 
the dashed line. 



the plasma parameters, indirectly coupling the thermophoretic effect to changes in the 
plasma parameters. 

5. Conclusion 

Using images of levitated dust particles, as well as of the plasma emission in the sheath, 
as well as electronic and Langmuir probe data, we have shown that thermophoresis has 
an effect on the discharge parameters in a capacitively coupled argon RF discharge. 
A slightly increased bulk resistance together with a more strongly reduced sheath 
reactance can be seen from the electronic data. Optical emission data seems to be 
consistent with this, in the sense that a clear change in the emission detected from the 
sheath region close to the electrode surface is observed, whereas the emission towards 
the plasma bulk does not change much with temperature. Langmuir probe data in the 
bulk is less conclusive, but seems to indicate a change in electron density, whereas the 
floating potential and electron temperature remain consistent with numerical estimates 
and do not show a clear trend with changes in the electrode temperature. 

Numerical modeling indicates that experiments running at constant pressure have 
a varying neutral atom density when thermophoresis is applied, due to an increase in 
the volume-averaged gas temperature. This is expected to have a noticable effect 
on the plasma parameters, consistent with our observations. Therefore, reports on 
thermophoresis experiments should clearly indicate if the DC bias was allowed to float, 
or was fixed, since the bias has an important effect on the dust levitation height and 
both types of experiments give different results for the equilibrium levitation height. 

In the past, such information was considered unimportant. Our data 
seems to suggest that possibly some aspects of experiments in complex plasma 
involving thermophoresis might have been overlooked. In some recent experiments, 
thermophoresis played an important role, but the effects discussed in this report were 
not explicitly considered. Discussions in these works involved, for instance, the correct 
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formulation of the thermophoretic force using dust particle levitation [S], the exact 
shape of the external confinement of particles levitated in a glass box with the help 
of thermophoresis [25 , as well as observations of oscillations in a single layer of dust 
particles at an unusually thin edge (due to additional thermophoresis) of a void, even 
though in that particular case the temperature gradient arose naturally [26]. Even 
though the relative changes in plasma parameters due to thermophoresis depend on 
the exact discharge settings, our results indicate that any explanation that assumes 
constant plasma parameters at different temperatures used to induce thermophoresis, 
are automatically invalid and can not describe the external confinement of dust 
particles in the discharge correctly for all temperatures. This is something that 
experimentalists in the field should be aware of. 
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